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The crossover between an impurity band (IB) and a valence band (VB) regime as a function 
of the magnetic impurity concentration in models for diluted magnetic semiconductors (DMS) is 
studied systematically by taking into consideration the Coulomb attraction between the carriers and 
the magnetic impurities. The density of states and the ferromagnetic transition temperature of a 
Spin-Fermion model applied to DMS are evaluated using Dynamical Mean-Field Theory (DMFT) 
and Monte Carlo (MC) calculations. It is shown that the addition of a square-well-like attractive 
potential can generate an IB at small enough Mn doping x for values of the p-d exchange J that are 
not strong enough to generate one by themselves. We observe that the IB merges with the VB when 
X ^ Xc where Xc is a function of J and the Coulomb attraction strength V . Using MC calculations, 
we demonstrate that the range of the Coulomb attraction plays an important role. While the on-site 
attraction, that has been used in previous numerical simulations, effectively renormalizes J for all 
values of x, an unphysical result, a nearest-neighbor range attraction renormalizes J only at very 
low dopings, i.e., until the bound holes wave functions start to overlap. Thus, our results indicate 
that the Coulomb attraction can be neglected to study Mn doped GaSb, GaAs, and GaP in the 
relevant doping regimes, but it should be included in the case of Mn doped GaN that is expected 
to be in the IB regime. 

PACS numbers: 71.10.-w, 75.50.Pp. 



I. INTRODUCTION 

The development of spintronics devices^ has moti- 
vated a large body of research on diluted magnetic 
semiconductors^ with the ultimate aim of creating ma- 
terials with Curie temperatures (Tc) above room tem- 
perature. This ambitious goal can only be achieved by 
a detailed understanding of the underlying mechanisms 
that govern the behavior of currently available DMS. 

Most theoretical approaches to study these materials 
start with one of two extreme regimes: (i) the limit of 
high Mn doping in which holes are directly doped into 
the valence band and, thus, are uniformly distributed in 
the sample (VB scenario and (ii) the limit of very 
low Mn doping in which holes are electrically bound to 
the impurity cores and an impurity band develops due to 
wave function overlap as the number of holes increases 
(IB scenario) ^ Researchers using the VB limit claim that 
it is valid for all the relevant dopings, namely x>l% 
in Gai_a;Mnj;As, and some experimental results support 
their viewi^i^ However, a similar claim is advanced by the 
groups promoting the IB scenario, i.e. that the IB exists 
up to the largest value of x that has been reached exper- 
imentally [x w 10%). This view also appears supported 
by the analysis of some experimental datai^i^ 

To solve this apparent puzzle, it is very important 
to study theoretically the DMS problem using unbiased 
techniques that provide reliable estimations for the value 
of X where the IB to VB crossover takes place. Such un- 
biased approaches could be provided by numerical tech- 



niques: in fact, the MC and DMFT methods have already 
been applied to a variety of phenomenological models for 
the DMS i^^d^'^^d^'^'^ These previous studies have been 
able to determine a crossover between the VB and IB 
behaviors, but only as a function of increasing values of 
the p-d exchange J. However, most experimental results 
appear to indicate that the realistic J for (Ga,Mn)As 
is approximately 1 e\ ^ which corresponds to the weak 
coupling regime in which no IB is generated by J alone. 
In fact, recent results obtained applying MC techniques 
to a six-orbital microscopic model, in which both the 
correct lattice geometry and the spin-orbit interactions 
were considered, indicate that (Ga,Mn)As is indeed in 
the VB regime for x > 3%^ In addition, DMFT tech- 
niques, which allow for the study of the very diluted 
(x <^ 1) regime, have shown that for values of J in the 
weak coupling region, an IB never develops as a function 
of a;, iiii2,i3,i4 jjowever, experiments based on electron 
paramagnetic resonance^ infrared spectroscopy^i^ and 
magnetization mcasurementsi^ of the electronic struc- 
ture of one Mn ion doped in GaAs have actually shown 
the existence of a shallow hole state with binding energy 
Eb=112.4 meV centered at the S=5/2 Mn ion. More- 
over, analytical studies indicated that Eb has contribu- 
tions from both the spin-dependent p-d hybridization and 
the Coulomb attraction between the hole and the Mn 
trapping center.— When additional Mn ions are added, 
the wave functions of the bounded holes will start to over- 
lap and an IB will develop. Further increasing x should 
widen the IB, locating it closer to the VB and eventually 
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a regime of complete hybridization with the holes doped 
into the VB is expected to occur. Thus, it is clear that a 
crossover from the IB to the VB regime should take place 
in (Ga,Mn)As as a function of x. 

In this paper, it will be argued that an IB-VB crossover 
will be missed in theoretical studies of materials with a 
weak J if the Coulomb attraction is disregarded, while 
materials with very strong J will be in the IB regime 
regardless of doping. In fact, here we explicitly show 
that by the simultaneous consideration of J and V in the 
formalism, the experimentally observed transition from 
IB to VB with increasing x can be understood. The or- 
ganization of the paper is the following: in Section |TT] 
the non-magnetic interactions in DMS are described; the 
model used and the DMFT technique are presented in 
Section Hm in Section HVl the results, including MC sim- 
ulations, are discussed, and Section |V] is devoted to the 
conclusions. 



II. SPIN-INDEPENDENT INTERACTIONS 
BETWEEN HOLES AND MAGNETIC 
IMPURITIES 

As remarked in the Introduction, most of the numerical 
work on DMS has been performed on models that focused 
on the role of the spin dependent p-d exchange J interac- 
tion between the spins of the localized impurities and the 
doped holes . This is certainly sufficient to capture 
qualitatively many of the properties of these compounds, 
including the generation of ferromagnetism. However, 
non-magnetic interactions between holes and impurities 
must be considered in order to improve the quantita- 
tive agreement with experiments. This additional po- 
tential term in the model has been generally referred to 
as "chemical disorder" (1/),— and it summarizes all the 
non-magnetic interactions between the localized impuri- 
ties and the holes. In this context, Tworzydlo^i used a 
short range potential (less than nearest-neighbors range) 
with a square- well form of depth Vb, and considered both 
positive (repulsive) and negative (attractive) values of Vq. 
The potential was introduced to explain an apparent x- 
dependence of the p-d exchange in Cdi_2;Mn2;S. Dietl^^ 
recently used the same approach to address apparently 
contradictory experimental results for Gai_a;Mn2;N. He 
also pointed out^ that this kind of extra potential term 
leds to a chemical shift in the standard impurity lan- 
guage, or to a valence-band offset in the alloy nomencla- 
ture, and that J and V are actually related^Si^i through 
the expression V/ J = 5(J7off-l-2erf)/4[/off where UcS is an 
effective correlation energy for the Sd shell, and is its 
energetic position with respect to the top of the valence 
band. However, the value of V is not easy to determine 
and, thus, it has been added as an extra free parameter 
by some authors (with V allowed to take both positive 
and negative values).— 12^>^ Other efforts focused just on 
the attractive Coulomb interaction between the holes and 
the impuritiesFi^i^i21 



Only some of the previously mentioned investigations 
have attempted to study the effects of the Coulomb at- 
traction at finite x with unbiased techniques. The au- 
thors of Ref. [2^ studied the case of a single Mn im- 
purity, considering the long-range Coulomb potential 
supplemented by a central cell correction with a gaus- 
sian or square-well shape, that is routinely introduced 
in calculations of bound state energies for impurities in 
semiconductors)^ For higher dopings, it is believed that 
the most important coulombic term is the central-cell 
contribution since the long-range potential is screened. 
In Ref. the coherent potential approximation (CPA), 
very similar in spirit to DMFT, was applied to a single 
orbital model which included both the spin dependent 
p-d hybridization J and an on-site central-cell Coulomb 
attraction V . It was claimed that the IB-VB crossover 
for (Ga,Mn)As using V~O.Q eV (chosen to reproduce, in 
combination with J=0.8 eV, the single impurity bound 
state energy) should occur for a; ~ 1 — 3%. In Ref. [IHl, a 
repulsive on-site potential was added. Both the repulsive 
and attractive cases were considered in Ref. How- 
ever, these important previous efforts did not present a 
systematic analysis of results as a function of J, V , and 
x, which is part of the goals of the present study. 

In this work we apply DMFT to a model that includes 
J and the Coulomb attraction V . The density of states 
(DOS) and Tq are studied in a wide range of couplings, 
hoppings, carrier fillings p, and Mn concentrations x, and 
estimations of the most appropriate values for different 
materials are made. We obtain the IB-VB crossover for 
a large class of DMS's and show that with a suitable 
strength V included, the IB regime can always be reached 
by decreasing the Mn concentration. 

III. MODEL AND DMFT FORMALISM 

The Spin-Fermion Hamiltonian used here and in sev- 
eral previous studies contains a kinetic f-term that de- 
scribes the hopping of holes between two neighboring i 
and j lattice sites {t is set to 1 to define the energy unit), 
an exchange interaction (EI) J^f-term that anti-aligns the 
carrier's spin with the magnetic moment of the impurity 
(considered classical) at site /, and a V-term that takes 
into account the on-site central-cell part of the attractive 
Coulomb potential,— 

^=-^ ^(cLc,.+H.c.)+2JH^SrS/-y^n,. (1) 

Here, {cia) is the creation (destruction) operator for 
a hole with spin a at site i, Si=c\^actj3Cip/2 is the hole's 
spin, S/=S'm/ is the classical spin of the local moment, 
and ni is the number of holes at /. 

Several details on the DMFT calculations were al- 
ready presented in Ref. [ll| for the case V=0, thus here 
only a brief summary is given and the modifications 
introduced by a non-zero V are remarked. DMFT uses 
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the momentum independence of the self-energy in infi- 
nite dimensions [E(p, ia;„)^E(ju;„), ujn=(2n + 1)'!tT^ 
and reproduces the physics of diluted correlated sys- 
tems in lower dimensionsi^ Within DMFT, the bare 
Green's function Go{ii-^n) contains all the information 
about the hopping of carriers onto and off magnetic 
(with probability x) and nonmagnetic (with proba- 
bility 1-x) sites. With ([1]) the full Green's function 
G{iujn) is solved by integration obtaining the result: 

{G{iuj^))^x{[G^\iuj^)^Jxnd+vi]-^) + {l-x){Go{iLo^)), 
where J=JhS.— This equation, complemented 
with the relation {Gq^ {iuJn))='Zn—{W'^ /16){G{iu!n)) 
valid within the assumption of a Bethe lattice^^ 
can be solved with a semicircular noninteracting 
DOS(a')=2ReV(W^/2)2-w2/7rW^ (z„=/^-f ia;„, fi is the 
chemical potential, and W=4:t is the bandwidth). Being 
spin diagonal, (Go) and (Gq ^) are expanded in powers 

of (Tz as: {a)=aoi+J2k '^k'^z , where ak^M'^, M being 
the order parameter used to detect the FM transition. 



To linear order in M we write {GQ^{iLLin))=B{iu!n)i 
+Q{iuJn)'^z and then Biiojn) is found from a 4-th order 
equation, 

W^[B± + V±JM] 1 
i?± = z„-x— p--^;j^-^-(l-x) — — , (2) 

that at /i=0 and with iw„^w gives us the low- 
temperature interacting DOS±{Lj)—~lm[B±{uj)\/-K for 
up (-f ) and down (— ) spin configurations.— The expres- 
sion for Q{iu!n)' 



Q+JM 



2J2Q/3 



16 i(B+F)2-j2 [(5+y)2_j2]2 



(3) 



leads us to an implicit equation for Tc in the form: 



4.tVK2J2B2 



[48B2 - 'i{\~x)W^]{[B+VY-.Py 



■ixW^B'^{\B+VY-J'^] - 2x^2 J2B2 
I 



= 1, 



(4) 



where B(iw„) is given by Eq. Q at M=Q. The Tc con- 
tained in LOn can be obtained from Eq. ^ numerically. 



IV. RESULTS 
A. General Analysis 

Let us start the discussion of results by considering the 
general dependence of a variety of quantities with the pa- 
rameters of the model. The DOS obtained from Eq. ([2|) 
at a;=0.035 is displayed in Fig. [T] for various values of J, 
M, and V. As observed in Fig. [TJa), the J-term alone is 
able to generate an IB but only if J /W exceeds a critical 
value Jc/VF~0.35. At reahstic couplings for (Ga,Mn)As 
(namely, J/W^=0.25 if we assume J«i~leV) there is no 
IB generated by the J-term alone. However, with the ad- 
dition of Coulomb attraction, when a value V/W^0.125 
is reached, then a well-defined split IB forms, as shown in 
FigHJb). No "symmetric" impurity band exists at high 
energies since the observed one is due to the carriers that 
are trapped in the vicinity of the core spins through the 
influence of V, and are fully aligned for M=l (Fig. [IJc)). 
The growth of J/W produces asymmetric low- and high- 
energy impurity bands if V^O (Fig. IHd)). 

We have observed that the coupling strength Jc/W for 
which the IB develops is a function of a;, namely the larger 
X is, the larger Jc/W becomes. Thus, wc used Eq. ([5]) 
to draw the phase diagram Jc/W vs. x at various val- 



ues of V. When V=Q the occurrence of an IB due only 
to the J-term requires a Jc/WkO.25 when x-^0, as seen 
in Fig. [2ta). When x^O and J/W<0.25 the addition 
of a potential V leads to the relation ( J-|-y)/W^«0.25 
to establish the boundary of the region where an IB de- 
velops. Our calculations also show that the boundary 
between the IB and VB regions in the full J-x plane just 
moves down by an amount A{V) after the introduction 
of the Coulomb attraction. This A{V) is independent 
of X indicating that Jc{x,V) = Jc{x,V=0) — A{V) as 
it can be seen in Fig. [2{a)j^ This means that an IB 
will be generated by a J < Jc(x, V=0) if a F such that 
(J+V) /W^Jcix) /W\v=o is added. Then, intuitively the 
effect of the addition of V is to renormalize J to a larger 
value. This result is not surprising because J has a dual 
effect: (i) it induces ferromagnetism, but (ii) it also tends 
to localize the holes near the impurity so that they take 
advantage of the antifcrromagnctic coupling. This last 
property is similar to the effect produced by the Coulomb 
attraction V. However, it would be expected that as x 
increases and more holes are added to the system, the 
wave functions of the holes will start to overlap, and 
as the holes become delocalized the effects of V should 
become less important. Thus, we would expect that the 
crossover boundaries between the IB and VB regions in- 
dicated in Fig. [IKa) should become closer to the V ^ 
curve as x increases, instead of remaining parallel as in 
the figure. Similar results have been observed in MC 
simulations 1^ We believe that the reason for this unex- 
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FIG. 1: (Color online) (a) DMFT low-temperature DOS at 
V=0, M=0, and different values of J/W. An IB forms if J/W 
exceeds a critical value «0.35. (b) DOS at M=0, J/W=0.25 
(believed to be realistic for (Ga,Mn)As), and different values 
of V/W. An IB forms if V/W^0.125. (c) same as in (b) but 
at V/W— 0.125 and for several values of M. The solid curve 
corresponds to DOS- while the dotted curve is for DOS+. 
(d) DOS at M=0, V/W=0.15, and various J/W. With a 
V/W^O the electron-hole symmetry is lost. In all frames the 
DOS is in arbitrary units and 1=0.035. At a;=0.05 we have 
reproduced the DOS obtained in Ref. [ij with CPA. 



pected behavior is related to the fact that here an on-site 
central-cell potential is being considered. This behavior 
can be corrected by considering a nearest-neighbor-range 
potentiai^^ or, within the DMFT framework, by consid- 
ering a phcnomcnological on-site potential that depends 
on X such as 



V{x) = VQexp{-{x/xof}, 



(5) 



where Xq can be roughly estimated using Mott's 
criterion^! as 



Xo 



0.25^ I ao 

as 



(6) 



with flo being the side of the cubic cell of the material 
and qb the Bohr radius for the bound impurity. For 
a material such as (Ga,Mn)As, which has an estimated 
as ~ SA, we obtain a;o=0.0014. The resulting boundary 
between the IB and VB regions is presented in Fig. [U^b) 
which indicates that for realistic values of J (0.2W) and 
Vo (O.IW) for (Ga,Mn)As, the crossover would occur for 
X < 0.5%. 

After having remarked that some paradoxes of the re- 
sults can be solved by extending the size-range of the 
attraction or, similarly, by reducing its strength with in- 
creasing X, here we will continue the discussion of the 
qualitative aspects for the case of the on-site central-cell 
potential. The main reason for it is to be able to com- 
pare our conclusions with previous results in the litera- 
ture since an on-site potential is the only approach used 



FIG. 2: (Color online) (a) The diagram J/W vs. x for vari- 
ous values of V. The solid curve defines the IB-VB crossover 
at V = 0. (b) The diagram J/W vs. x for (Ga,Mn)As 
when V is X dependent. The Vj^O curves all join at a::~0.005, 
that marks the Mn doping concentration beyond which the 
Coulomb attraction is no longer relevant and the IB disap- 
pears for realistic couplings, (c) The diagram V/W vs. x 
at a realistic ratio for (Ga,Mn)As J/W=0.25 with an on-site 
Coulomb attraction. The inset shows the merging of the im- 
purity and valence bands with increasing x, at V/W=0.066. 
(d) The diagram J/W vs. V/W at various x. The inset 
shows the DOS at J/W=0.2, V/W=0.U8, and a;=0.03. Since 
Jc/W is a;-dependent, the VB "triangle" shrinks (expands) 
as X decreases (increases), with the shrinking saturating at 
Jc/W ^0.25. 



in previous numerical investigations! There are still 
some quantitative aspects that may need the x depen- 
dent potential of the previous paragraphs, and those will 
be clarified below. 



Focusing on the on-site potential, it can be observed 
that even if J/W<Jc/W, the IB regime can in general 
be reached either by increasing V at fixed x, or by de- 
creasing X at fixed V (see Fig. [5Jc)). While at x^O the 
carriers trapped due to V in the vicinity of each Mn core 
spin reside in an impurity-like bound state, as x increases 
the wave functions that describe the bound state at the 
manganeses start overlapping (due to the combined ef- 
fects of V and J) producing an IB that at a critical Xc 
merges with the VB. The renormalization condition ob- 
tained in our calculations yields an IB-VB boundary in 
the diagram J/W vs. V/W, for a fixed X , clS shown in 
Fig. [DJd). This boundary deviates from linear only for 
very small values of J/W which is not a physically in- 
teresting region. According to the results in Fig. [Hd) 
the area of the VB region is a minimum for a; ^ and 
increases with increasing x. 
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B. Specific Results for (Ga,Mn)As and Otlier 
Compounds 

The literature does not provide a unique value of V for 
the case of (Ga,Mn)As. The main reason is that the value 
of V necessary to generate a bound state upon doping by 
one hole is a function of both J and the bandwidth W, as 
it can be observed from the results presented in Table [J 
Thus, in Ref. [l^l a value of V =2.3 eV is determined for 
J=0.9 eV with W 10 eV since a Luttinger-Kohn en- 
ergy band is used, while in Ref. [isl. V=0.6 eV is used with 
J— 0.8 eV and W=A eV. In both cases, V is determined 
by requesting that for a single impurity doping a bound 
state at i?b= 112 meV appears as the combined result of 
the magnetic and Coulomb interactions. Our calculations 
indicate that the parameters of Ref. [1^ provide an IB- 
VB crossover at 0.5% while we recovered the value 
Xc ~ 3% of Ref. [1J| using the parameters that they pro- 
vided. The discrepancy shows that the values assumed 
for W and J play an important role in the determination 
of V and Xc- The expression given by Dietl?^ provides an 
estimation of the non-magnetic impurity potential that 
may include more than Coulomb interactions. It is eval- 
uated using experimental data. For x ~ 7 %38,39 .^j^^^ 
W = 3eV and J=l eV, the ratio |y/ J|=0.55 is obtained. 



The potential turns out to be repulsive y=— 0.55 eV. No- 
tice that while the estimations of V performed for x — > 
provides positive values, compatible with an attractive 
potential, the estimations at finite doping do not. As 
pointed out in the previous section, this indicates that it 
may be necessary to use an x-dependent expression for 
the non-magnetic interactions. 

The phenomcnological potential proposed in Eq. ([5]) 
will provide an IB-VB crossover at a::^0.1% for all the at- 
tractive values of V provided above, as seen in Fig. [21^b). 

We can make estimations of Xc for (Ga,Mn)As and for 
other Mn doped III-V materials as well. The value of J is 
expected to be inversely proportional to the volume of the 
cubic cell of the material Qq, according to the chemical 
trends, and the energy of the bound state for one Mn 
impurity has been measuredi^ From these data, we can 
estimate V for different values of W, with results given 
in Table IH that also includes ao for each material and 
the estimated value of ae ~ h/\/2'mkEi, where rrik ~ 
'me/ ill — (673 -l- 472)75) with TTie the electron mass and 
7i the Luttinger parameters4S Then xq can be obtained 
from Eq. ^ and is also shown in the Table. Xc (ic) 
indicates the estimated values of the doping for which 
the IB-VB crossover occurs for an on-site (x-dependent) 
potential {V{x) given by Eq. ([5])). 



TABLE I: DMFT calculated values of V that produce a bound state with energy Ei, for the values of J and bandwidth W 
shown corresponding to the indicated DMSs. The calculated doping density Xc (ic) at which the IB/VB crossover occurs for 
an a;— independent (dependent) potential is listed. The IB label indicates that the material is in the IB regime at all x£{0, 1]. 
Values of ao, as, and xo (see text) for each material are also shown. 



Material 


J (eV) 


Et (eV) 


W{eV) 


V(eV) 


Xc{%) 


ao (A) 


as (A) 


Xo 


Xc (%) 








10 


2.7 


IB 








7.2 


(Ga,Mn)N 


2.5 


1.4 


8 
6 
4 


2.014 
1.31 
0.47 


IB 
IB 
IB 


4.42 


1.6 


0.082 


9.3 
21 
IB 








10 


2.4 


5.2 








0.422 


(Ga,Mn)P 


1.34 


0.41 


8 
6 
4 


1.786 
1.173 
0.525 


8.3 
16.7 
30 


5.45 


4.5 


0.007 


0.493 
0.637 
2.14 








10 


1.883 


0.52 








0.059 


(Ga,Mn)As 


1.2 


0.112 


8 
6 
4 


1.324 
0.761 
0.19 


0.85 
1.35 
3.1 


5.65 


8 


0.0014 


0.068 
0.09 
0.37 








10 


1.74 


0.025 








0.00044 


(Ga,Mn)Sb 


0.96 


0.016 


8 
6 
4 


1.232 
0.698 
0.175 


0.045 
0.064 
0.13 


6.10 


39 


0.00015 


0.00053 
0.00065 
0.0014 



It is clear that for all relevant values of x, (Ga;Mn)As 
is in the VB regime. The crossover, for realistic values 
of W, occurs at x<l% for both on site and cc-dependent 
potentials. Thus, even including the Coulomb attraction. 



our results indicate that the IB regime is not expected 
to play a relevant role in this material. A similar picture 
emerges for (Ga,Mn)Sb. In this case the IB-VB crossover 
is expected to occur for such small values of impurity dop- 
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FIG. 3: (Color online) (a) Tc vs. p at V/W = for several 
values of J/W . (b) Tc vs. p at J/W=0.25 for various values 
of V/W. In both frames a;=0.05. 



ing that for all practical purposes the Coulomb attraction 
can be neglected. 

On the other hand, the IB regime seems to dominate 
the physics of (Ga,Mn)N. Considering J— 2.5 eV, within 
our model we found that even for the largest value of W 
considered (namely, VF=10 eV) J/W is strong enough to 
generate an IB region below some finite Xc(W), even if no 
Coulomb attraction is considered. However, since the sin- 
gle hole bound energy for GaN is 1.4 eV, i.e. much larger 
than the 0.113 eV value observed in GaAs, it is clear that 
the Coulomb-attraction term has to be incorporated. In 
the table we show the values of V that together with J 
will produce the bound state for different values of the 
bandwidth W. Our calculations show that with an on- 
site potential (Ga,Mn)N will be in the IB regime for all 
relevant values of x (we studied up to x = 80%). This 
is still true when an x-dependent V is considered since 
even in the case for the largest bandwidth considered the 
crossover is expected to occur at x ~ 7.2%. Coulomb 
attraction should therefore be included to study this ma- 
terial. 

Our results for (Ga,Mn)P indicate that despite the 
deeper position of the bound state in the gap, studies 
neglecting the Coulomb attraction could be performed, 
particularly for x > 3%. 

For completeness, and to compare with previous 
calculations we present the Tc vs. p dependence 
obtained from Eq. ^ at a-=:0.05, for different values 
of J's and no Coulomb attraction in Fig. El^a). For 
J/W<^Jc/W, Tc is low and almost independent of p. 
When J/W>Jc/W, i.e. in the IB regime, Tc vs. p 
is semicircular with a maximum at p=x/2, in agree- 
ment with previous results for one-orbital models The 
behavior of Tc vs. p at different values of V/W for 
J/W = 0.25 is shown in Fig. [3]Jb). Comparing with 
the curves in part (a) of the figure it is clear that V 
increases the effective value of J. Our results agree with 




FIG. 4: (Color online) (a) Tc vs. J/W at p=0.015 and x = 
0.05 calculated with DMFT for different values of V/W. (b) 
Tc vs. J for different values of V at ph=0.3 and x — 0.25 
obtained by MC. (c) The density of states (DOS) for J/t = 1 
and V = (black line); for an on-site Coulomb attraction 
V — 2 (red line; the curve has been shifted vertically for 
clarity); and for a finite-range Coulomb attraction with on 
site intensity V and next nearest neighbors intensity V'=V=2 
(dashed green line). The vertical lines indicate the chemical 
potential. For clarity, the curves for finite Coulomb attraction 
strength have been shifted along uj so that the central peak 
in the DOS of all the curves coincides. 



Ref. |l3| and confirm that an on-site square- well V simply 
renormalizes J. The dependence of Tc on J for different 
values of V is shown in Fig. Ufa). V boosts Tc at small 
and intermediate J/W, while at large J/W's no change 
is observed because within DMFT the Tc saturates as 
J oo. However, as it will be discussed in the following 
section, we believe that the renormalization of J for the 
physically relevant values of x, such as the one used in our 
figures, is an artifact of the on-site range of the Coulomb 
attraction and, thus, we do not expect it to play a role 
in enhancing the Tc of real materials. 



C. Monte Carlo Simulations 

Hamiltonian ([1]) was also studied here using a real- 
space MC technique with the Mn core spins treated clas- 
sically. Details are not provided since the technique has 
been widely discussed before in the context of studies of 
manganites4i The simulations were performed using cu- 
bic lattices with 4"^ sites at a::=0.25. Finite-size effects 
have been monitored by running some points on 5'^ clus- 
ters. A random starting spin configuration has been se- 
lected as the starting point for each temperature T. The 
spins were allowed to evolve for a total of 10^ MC steps, 
with the first 5x10"* steps being discarded to thcrmalize 
the starting configuration. 

At J/t=l, and ph=p/x=0.3, a value V^l for the on- 
site Coulomb attraction increases Tc by as much as 33%, 
as shown in Fig. IHJb). This agrees qualitatively with 
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the DMFT results. The figure shows clearly how V ef- 
fectively " renormalizes" J. Since the curve Tq vs. J for 
V=0 has a maximum at J"*"^ the effect of V is to increase 
Tc for values of ( J + T/)< J"""^, while Tc decreases with 
V for values of (J + V) > J™"^. Although this renormal- 
ization has been previously reported,— we do not believe 
that it will play a role in the relevant range of doping for 
most DMS. As wc pointed out in subsection IIV Al the 
on-site range of the Coulomb attraction induces unphys- 
ical behavior by exaggerating hole localization for values 
of X for which overlap of the hole wave functions should 
occur. While finite range attraction cannot be studied 
with DMFT, it can be done with MC simulations but 
at the price of not being able to access the low doping 
regime at which the IB-VB crossover would be expected 
to occur for a material such as (Ga,Mn)As. 

In Fig. HJc), we present the DOS obtained with MC 
for J/t=l and V=Q for x=25% indicated by the black 
continuous line. The peaks are due to the finite size of 
the system, and each of them can be identified with the 
spikes that appear in the DOS of a non-interacting sys- 
tem in the same lattice. Thus, at this value of J, there 
is only a VB in the DOS, i.e., the magnetic interaction 
is not strong enough to develop an impurity band. The 
position of the chemical potential fj, is indicated by the 
black dashed line. Upon adding an on-site Coulomb at- 
traction V=2, we observe that an IB develops as indi- 
cated by the red line in the figure, that has been shifted 
upwards along the vertical axis for clarity. This IB is due 
to the localization of the holes induced by the on-site po- 
tential. The chemical potential denoted by the dotted 
red line indicates that only states in the IB are occupied. 
However, when the range of the potential is increased to 
next-nearest neighbors, as indicated by the green dashed 
line in the figure, it can be seen that the IB dissapears 
although the intensity of the potential has not changed. 
This occurs because, at this large doping, the extended 
potential allows for a more uniform distribution of the 
holes. As it can be seen in the figure, the DOS for V~Q 
and for finite extended V have an almost perfect overlap. 
This shows that the use of on-site Coulomb attraction po- 
tential can lead to missleading results and authors have 
to be cautious when using this approximation. 

V. CONCLUSIONS 

Our combined DMFT-MC study shows that the 
Coulomb attraction by acceptors needs to be considered 
to obtain correctly the IB-VB crossover as a function 
of impurity doping concentration x in models for DMS. 
However, for most materials we find that the crossover 



occurs at very low levels of doping, outside the regime 
in which high Tq would be expected. We also find that 
a doping-independent on-site square-well potential acts 
as a renormalization of the coupling J in an extended 
doping range up to x = 80%. However, this appar- 
ent boost to the J-term at all Mn dopings is unphysi- 
cal, since the effect of V should be x-dependent beyond 
some critical value. Our MC simulations demonstrate 
that this a;-dependence is achieved naturally by consid- 
ering a longer range (next-nearest neighbors) square- well 
attraction, which is beyond the capability of the single 
site DMFT which can deal with on-site interactions only. 
Thus, a phenomenological a;-dependent Coulomb attrac- 
tion was introduced. With this modification, we have 
shown that for (Ga,Mn)As, the Coulombic attraction V 
infiuences the physics of the material only at small Mn 
doping , i.e a;<0.5%. This result shows that it is cor- 
rect to apply theories that consider the J-term only for 
studying the properties, including the Curie temperature, 
of (Ga,Mn)As at the relevant values of Mn concentra- 
tions a;~l%-10%. On the other hand, we found that the 
Coulomb attraction will play a relevant role, and should 
be included, in studies of Mn-doped GaN. 

Summarizing, here we have shown that the addition 
of an attractive Coulomb potential is the necessary in- 
gredient to explain the transition from the IB to the VB 
regime as a function of Mn-impurity doping concentra- 
tion in materials for which the magnetic interaction J is 
not strong enough to bind a hole. However, we find that, 
except for the case of (Ga,Mn)N, the crossover occurs 
at very low doping in a regime in which high ferromag- 
netic critical temperatures would not be expected and, 
thus, the effective value of J will not be affected. As a 
consequence, it is not necessary to include the Coulomb 
attraction in the calculations. In addition, we show that 
an on-site attractive potential does not capture the over- 
lap of localized hole wave-functions that should occur as 
a function of doping and it provides unphysical results. 
Thus, to study materials such as (Ga,Mn)N, in which the 
Coulomb attraction is relevant, a nearest-neighbor finite 
range potential has to be used. 
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